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A standard quarter-wavelength multilayer optical coating will produce the highest reflectivity 
for a given number of coating layers, but in general it will not yield the lowest thermal noise 
for a prescribed reflectivity. Coatings with the layer thicknesses optimized to minimize thermal 
noise could be useful in future generation interferometric gravitational wave detectors where coating 
thermal noise is expected to limit the sensitivity of the instrument. We present the results of 
direct measurements of the thermal noise of a standard quarter-wavelength coating and a low noise 
optimized coating. The measurements indicate a reduction in thermal noise in line with modeling 
predictions. 

PACS numbers: 



INTRODUCTION 

Several second generation gravitational wave detectors 
are under construction now (Advanced LIGO 1, 2], Ad- 
vanced Virgo, GEO-HF Q). In general, the sensitivity of 
these instruments will be limited at low and high frequen- 
cies by quantum noise, however, in a critical mid-band 
around 100 Hz, the thermal noise of the test mass mirror 
coatings is a significant limit. A typical sensitivity curve 
for a second generation detector is shown in Fig. Q~] In 
light of this fact, the study of thermal noise in dielectric 
coatings is an active area of research in the gravitational 
wave community. 

The optical cavities at the heart of the length-sensing 
mechanism of gravitational wave interferometers use mir- 
rors made with multilayer dielectric coatings to produce 
the high reflectivities that are required. For good per- 
formance (low absorption) at Ao = 1064 nm (the operat- 
ing wavelength of most gravitational wave detectors), the 
coatings are usually made of alternating layers of silica 
(Si02) and tantala (Ta205). In an elegantly conceived 
experiment that studied mechanical loss in these coat- 
ings, Penn et al. Q showed that the primary source 
of mechanical dissipation in these coatings was in the 
tantala layers rather than in the silica layers or at the 
interfaces between layers. This result suggests two possi- 
ble mechanisms for reducing the mechanical dissipation 
and, therefore, the in-band Brownian noise |5|. First, one 
can alter the chemistry of the tantala to reduce its me- 
chanical loss. This has been accomplished by doping the 
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FIG. 1: Projected noise floor of Advanced LIGO Q. As it 
stands, Brownian noise in the test mass coatings will prevent 
the instrument from reaching the quantum limit. 



tantala with titania (TiCh) to good effect [6j. Second, 
one can modify the geometry of the coating to reduce 
the total amount of tantala while preserving the reflec- 
tivity of the coating. Such a coating is referred to here as 
an "optimized coating" and is the subject of this paper. 



A standard high-reflectivity multilayer dielectric coat- 
ing consists of alternating layers of high and low index of 
refraction materials where the thickness of each layer is 
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1 /4 the local wavelength of the light. This design requires 
the minimum number of layers to achieve a prescribed re- 
flectance. It does not, however, yield the lowest thermal 
noise for a prescribed reflectance. If there is a difference 
in the mechanical loss of the high and low index mate- 
rials then the overall coating dissipation can be reduced 
by decreasing the total amount of the high loss material. 
We have developed a systematic procedure for designing 
minimal-noise coatings featuring a prescribed reflectivity. 
We have manufactured mirrors based on such a design at 
the Laboratoire des Materiaux Avances and measured 
the broadband noise floor of the Thermal Noise Interfer- 
ometer (TNI) at the California Institute of Technology 
with these mirrors installed. We then compared this to 
an earlier measurement of the noise floor of the TNI when 
mirrors with standard quarter-wavelength (QWL) coat- 
ings manufactured at Research Electro-Optics, Inc. were 
in place. This paper summarizes the relevant theoretical 
background, and presents the results of these measure- 
ments. 



THE THERMAL NOISE INTERFEROMETER 

Until 2002 Q , surface fluctuations of thermal origin in 
the mirrors of an optical cavity had never been directly 
observed. Up to that point the small scale of the fluc- 
tuations meant that they were always below the level of 
other noise sources such as laser frequency noise, shot 
noise, etc. The TNI is a test bed interferometer specif- 
ically designed to detect these surface fluctuations. It 
consists of a Lightwave Model 126 NPRO (Non-Planar 
Ring Oscillator) laser, a triangular mode cleaner cavity 
used to spatially filter and stabilize the frequency of the 
laser beam, and two high finesse test cavities where cav- 
ity length noise measurements are made. The key fea- 
tures of the TNI design that enable resolution of thermal 
noise are a relatively high power laser of nearly 1 W, fre- 
quency stabilization by the mode cleaner by a factor of 
1000 compared to the free-running laser frequency noise 
in the detection band, a relatively small beam radius of 
164 /im to amplify the effect of the surface fluctuations, 
short test cavities to reduce the laser frequency stabi- 
lization requirement, and two identical test cavities to 
permit common mode noise rejection. See reference Q 
for more details about the TNI. 



The layout of the TNI is shown in Fig. d All three 
optical cavities are under vacuum and every cavity mir- 
ror is suspended like a pendulum from a loop of steel 
wire to preserve its high mechanical Q and to isolate it 
from external sources of noise. The Pound-Drever-Hall 
technique 13, II | is used to keep the cavities locked on 
resonance. In the case of a test cavity, a signal is gener- 
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FIG. 2: Schematic of the Thermal Noise Interferometer. 
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FIG. 3: Block diagram of the feedback loop of a test cav- 
ity. I and 5v denote fluctuations in cavity length and laser 
frequency, respectively. Data is taken at the point V . 



ated proportional to the deviation from the cavity reso- 
nance and this signal is fed back, after suitable amplifi- 
cation and filtering, to electromagnetic actuators at the 
end mirror of the cavity to control the cavity length. A 
block diagram of a test cavity servo is shown in Fig. [3] It 
is a simple matter to show that a fluctuation measured 
at the readout point, V , can be converted to its equiv- 
alent fluctuation in cavity length, £, with the following 
formula, in the spectral (frequency) domain: 



T~ 1 V 



1 - DCMH 
DC 



V. 



(1) 



It is imperative that the (spectral) transfer functions of 
the servo elements D, C , M, and H are known accurately 
so that this conversion is valid. This calibration error is 
the largest source of uncertainty at the TNI. 



In a single test cavity at the TNI, laser frequency noise 
begins to dominate the noise floor above about 6 kHz. 
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FIG. 4: Block diagram illustrating how the difference signal 
readout, Vd, is related to the length noise in the "North," In, 
and "South," Is, test cavities. An and As are attenuators 
used to balance out any difference in the responses of the 
servos Tjv and Ts- 



The frequency range where thermal noise is dominant 
can be greatly extended if common mode rejection is im- 
plemented between the two test cavities. This is done 
by using a commercial pre-amplifier (Stanford Research 
Systems SR560) to subtract the readout signal of one test 
cavity from the other as shown in Fig. |4j If the two test 
cavity servos are well matched, the difference between the 
readout signals, Vd, is related to the difference between 
the cavity length fluctuations in the following manner: 



e d = T^V d = - (A N ^ + A-/T s ') V d . 



(2) 



We refer to the two test cavities as the "North" and 
"South" cavities and the subscripts N and S reflect this. 
T^ 1 and T<7 1 are the transfer functions given by Eq. (1) 
for the North and South cavities, respectively. By con- 
verting the difference readout in this manner, we can 
compare the noise spectrum of the TNI to the theoreti- 
cally predicted thermal length noise. 



Once the common mode noise is removed from the TNI 
by subtracting the test cavity readouts from each other, 
the noise floor of the instrument above about 20 kHz be- 
comes dominated by shot noise. For each test cavity, the 
magnitude of the shot noise is determined by the power 
incident on the cavity's photodetector. This is light that 
is reflected from the test cavity and it consists primar- 
ily of the component of the test cavity input beam that 
is not mode-matched to the cavity (the mode-matched 
component is transmitted through the cavity when it is 
resonant). The power depends on how well the cavities 
are aligned with their input beams and this can vary be- 
tween measurements. The shot noise as a function of 
power on the photodetectors has been measured using 
a heat lamp as a shot noise limited source of radiation. 
Accordingly, determining the magnitude of shot noise at 
the time a thermal noise measurement is made is sim- 
ply a matter of noting the power on the photodetectors 
at that time and using the results from the heat lamp 
measurement. 



OPTIMIZED COATING DESIGN 

Designing an optimized coating is a matter of finding 
the configuration of layer thicknesses that minimizes the 
Brownian noise of the coating as a whole for a prescribed 
reflectance. Since the QWL design has the maximum 
reflectance for a given number of layers, the reflectance 
can only be maintained by adding more layers to the 
coating. To design the optimized coating, a means of 
calculating the reflectance of a coating as a function of 
the layer thicknesses is needed. 

For time- harmonic (e J ) normal plane wave incidence, 
the electric and magnetic fields at the i-th interface of a 
coating are related to the fields at the (i + l)-th one by 
the propagation matrix of the i-th layer, Mj [lij ]. 



At the TNI, thermal noise in the coatings can be ob- 
served for over a decade of frequencies starting at around 
800 Hz. This is above the unity gain frequency of the 
servos (typically around 500 Hz) where the conversion 
formula simplifies to: 
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(3) 



Here we have used the fact that the servo element C, 
which is the ratio of the frequency of the laser to the 
length of the cavity, is the same for both test cavities. 
All the servo responses in Eq. (3) are frequency inde- 
pendent constants. The elements An and As are at- 
tenuators whose gains are adjusted so that the product 
AnDn = A$Ds, to maximize the common mode rejec- 
tion ratio. 
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where Zq = y //q /eq is the characteristic impedance of 
the vacuum, m is the refractive index of the layer, and 
5i = kili is the phase thickness of the layer. It is conve- 
nient to write the phase thickness as Si = 2irzi, where Zi is 
the thickness of the layer in units of the local wavelength, 
Zi = k/Xi = lifii/XQ. See Fig. [5] for an illustration of a 
multilayer coating. The input impedance of the coating, 
defined as the ratio of the electric and magnetic fields at 
the vacuum-coating interface, Z- m = E\/H\, can then be 
obtained by chain multiplying the propagation matrices 
of all the layers in the coating: 
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FIG. 5: Diagram of a M layer coating. U is the thickness 
of the i-th layer and n; is its index of refraction. Ei and Hi 
are the electric and magnetic fields, respectively, at the i-th 
interface. 
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In the last step we used the fact that the fields at the 
coating-substrate boundary can both be written in terms 
of the electric field of the transmitted beam, E t . The 
reflection coefficient of the coating, T, in terms of the 
input impedance is: 
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(6) 



The (power) reflectance of the coating is R = \T\ 2 . With 
this, we have a means of calculating the reflectance of a 
coating from the layer thicknesses. We also need a means 
of calculating the Brownian noise of the coating from the 
layer thicknesses. 



The power spectral density (PSD) of the Brownian 
noise of a mirror is given by the following formula fl3|: 
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where ks is Boltzmann's constant, T is the temperature, 
w is the half-width of the Gaussian laser beam, a s is 
the Poisson ratio, and Y s the Young's modulus of the 
substrate. The PSD is proportional to the effective loss 
angle of the mirror, (p e s, and, if the substrate has a high Q 
factor, this is dominated by the loss angle of the coating, 
so that o ff <j) c . A precise formula for the coating loss 
angle is very complex (l3| . nonetheless, in the limit of 
small Poisson's ratios, it can be written as a weighted 
sum of the thicknesses of the materials in the coating: 



= b L d L + b H d H , 



(8) 



where dt, t B is the total thickness (in units of local wave- 
length) of all the low (L) and high (H) index layers, 



which, in our case, are made of silica and tantala, re- 
spectively. The weighting factors are given by: 
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where, Y and (f> denote the Young's modulus and loss 
angle and the subscripts L, H, s refer to silica, tantala, 
and substrate, respectively. Unfortunately, there is much 
uncertainty in the values of the loss angles for thin-film 
materials so the ratio 7 = bujb^ can only be said, with 
confidence, to lie somewhere between 5 and 10. The op- 
timized design was chosen among a few alternative ones 
yielding close-to-minimum noise, as the least sensitive to 
the inherent uncertainty in the value of 7 as regards the 
noise PSD reduction. 

The problem of finding the optimal configuration of 
coating layer thicknesses may seem a difficult one if the 
thicknesses of each and every layer were allowed to be 
free variables. However, preliminary simulations based 
on genetic minimization of the coating loss angle at pre- 
scribed reflectance, using two refractive materials and 
no prior assumptions about coating geometry, support 
the conclusion that the optimal configuration is periodic 
and consists of a stack of identical high/low index layer 
pairs or "doublets," with the exception of the terminal 
(top/bottom) layers [3]. Accordingly, in designing our 
optimized coatings we restrict to only these structurally 
simple stacked doublet configurations with tweaked ter- 
minal layers. 



Further simplification is possible, though. Figure|5]is a 
plot of contours of constant reflectance and constant ther- 
mal noise (in the case where 7 = 7) as functions of the 
thickness (in units of local wavelength) of the silica (zl) 
and tantala (zh) layers for a coating consisting of four 
identical tantala/silica doublets (for a total of eight lay- 
ers). In the figure, the reflectance of the coating increases 
to a maximum at the point (zl,Zh) — (0.25,0.25), the 
QWL configuration. The figure indicates that, for a 
given reflectance, there is not much difference in ther- 
mal noise between the point of minimum thermal noise 
(the white dot in the figure) and the point along the same 
reflectance contour where zl + zh = 0.5 (the black dot). 
Actually, as the number of doublets is increased, the re- 
flectance contours become even more squeezed along the 
half-wavelength doublet dashed line and the difference 
in thermal noise between the above two points becomes 
negligible. Thus, to simplify matters further, we focus on 
configurations where the doublets are a half-wavelength 
thick. This leaves only two free parameters to adjust: 
the number of doublets, N, and the thickness of the tan- 
tala layers, zh, the thickness of the silica layers being 
zl = 0.5 — zh- The coating design then proceeds as 
follows: 
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FIG. 6: Contours of constant reflectance (ellipse-like curves) 
and thermal noise (straight lines) for a stack of four identi- 
cal tantala/silica doublets as functions of the silica (zl) and 
tantala (zh) layer thicknesses (in units of local wavelength). 
The arrow pointing from darker to lighter areas indicates the 
direction of increasing thermal noise. In this plot 7 = 7. The 
dashed line consists of points where zl + zh = 0.5. The 
white dot marks the point of minimum thermal noise on that 
particular reflectance contour. The black dot marks a point 
on the same reflectance contour that satisfies the condition 
zl + zh = 0.5. This point has slightly greater thermal noise 
than the minimum noise point but the difference becomes 
negligible as the number of doublets (4 in this case) increases. 



1. Starting with the standard QWL design with 
Aqwl doublets and reflectance -Rqwl, use equa- 
tion (8) to calculate the loss angle of the coating. 

2. Add one doublet to the coating then reduce the 
thickness of the tantala layers (since they are more 
lossy) and increase the thickness of the silica lay- 
ers (to maintain the half-wavelength doublet con- 
dition) until the desired reflectance, Rqwl (calcu- 
lated from Eqs. (4)-(6)), is recovered. Recalculate 
the loss angle of the coating. 

3. Repeat step 2 and keep repeating until a design 



with a minimum coating loss angle, 
tained. 
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FIG. 7: The structure of the standard quarter wavelength 
coating. 
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TABLE I: The material parameters used in the design of the 
optimized coating. 



coating is Rqwl — 0.9997225 (the transmittance is 278 
ppm). Following the procedure outlined above, we de- 
signed an optimized coating with the same reflectance. 
The material parameters used in the design process are 
shown in Table |TJ Figure [8] is a plot that shows how the 
coating loss angle varies as the number of doublets is in- 
creased. As the figure indicates, the minimum loss angle 
design consists of 17 identical doublets. As suggested by 
the genetic simulations, this design was improved upon 
by tweaking the thicknesses of the end layers to further 
reduce the thermal noise while keeping the reflectance un- 
changed. The final optimized design is shown in Fig. |9] 
If 7 lies somewhere between 5 and 10 then the ratio of 
the coating loss angles, r = ^opt /</>qwl , should lie some- 
where between 0.876 and 0.817. Assuming 7 = 7 (the 
value used in the design process), then r — 0.843. 



The reference QWL coating measured with the TNI 
consists of 14 silica/tantala doublets and is shown in 
Fig. [7J For protective purposes, the topmost layer of the 
coating is silica for its hardness. It is half-wavelength, 
so as to have no influence on the coating transmission 
properties, as seen from Eq. (4). The reflectance of the 
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FIG. 8: Plot of the coating loss angle (normalized by the 
loss angle of the QWL coating) as a function of the number 
of doublets, N. The corresponding values of the ratio zh/zl 
are shown at the top of the plot. The QWL and the minimum 
noise (optimized) coatings are indicated in the plot. 
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Frequency (Hz) 

FIG. 10: Plot of the spectral density of cavity length noise for 
the QWL coating along with the shot noise (dashed line) and 
the Brownian noise (solid line with f~ 1 ^ 2 frequency depen- 
dence). To extract the loss angle of the coating the incoher- 
ent sum of the Brownian and shot noise (solid black curve) is 
fit to the data between the vertical lines at 2 kHz and 4 kHz 
with the loss angle as the free parameter in the fit. 

kHz the shot noise begins to dominate the noise spec- 
trum. The large peaks above 20 kHz are body modes of 
the mirror substrates. To extract the loss angles of the 
coatings, we focused on the region around 3 kHz since it 
is approximately at the center of the coating noise dom- 
inated region. 
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FIG. 9: This is the structure of the optimized coating. This 
coating was designed for minimal Brownian noise and a trans- 
mittance of 278 ppm. 



MEASUREMENT RESULTS 

To obtain the length noise in the test cavities at the 
TNI, we used a spectrum analyzer (Stanford Research 
Systems SR780) to measure the power spectral density 
(PSD) at the instrument readout, Vd- After calibrating 
the instrument by measuring the responses of the servo 
elements in Fig.[3J we converted this readout PSD (actu- 
ally the square root of the PSD) to its equivalent length 
noise spectrum using Eq. (2). The resulting length noise 
spectrum of the QWL coating is shown in Fig. [TO] The 
spectrum of the optimized coating is similar. Note that 
from about 500 Hz to 20 kHz the measured length noise 
has a slope of f^ 1 ^ 2 that is characteristic of Brownian 
noise (the square root of the PSD in Eq. (7)). Above 20 



In the process of taking measurements of the noise floor 
of the instrument, we observed a certain degree of insta- 
bility that is mainly due to instability of the optical gains 
of the servos, represented by the element D in the servo 
diagram in Fig. [3] To mitigate this, we took multiple 
measurements for each coating so that we could average 
the results. For each coating, we measured eight spectra 
from 2 to 4 kHz. Of these eight, four had a frequency bin 
width of 4 Hz, for a total of 500 points in each spectrum, 
and four had a 16 Hz bin width, for a total of 125 points 
in each spectrum. For each spectrum, we performed a 
least-squares fit of the function: 

i(f) = V^s B (f) + s s (io) 

where is the PSD of Brownian noise (Eq. (7)), and 

S s is the PSD of shot noise, which is obtained from the 
heat lamp measurement as described earlier. The coating 
loss angle is the only free parameter in the fit. The factor 
of 4 in the equation above is due to the fact that there 
are a total of four mirrors in the two test cavities, so the 
total PSD of Brownian noise is four times that of a single 
mirror. The resulting loss angles are shown in Table |H] 
and Fig. [TT]is a graphical representation of these values. 

The mean loss angles are </>q-wl = (8.43 ±0.35) x 10~ 6 
for the QWL coating and t = (6.90 ±0.18) x 10~ 6 for 
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TABLE II: Table of the coating loss angles (xlCT 6 ) that give 
the best fit to the data along with the mean value. *Mean 
after discarding anomalously high values. 



the optimized coating. To be conservative, we used the 
standard deviations of the two sets of eight loss angles for 
the errors of their respective means rather than the stan- 
dard deviations divided by the square-root of eight (the 
number of measurements for each coating). This was 
done because at least some of the variation in the loss 
angles of each coating arises from the instability in the 
optical gain mentioned earlier and this variation is non- 
Gaussian. The ratio of the means is r — 0.82 ± 0.04. Re- 
call that the predicted ratio of optimized to QWL coating 
loss angles (assuming 7 = 7) was 0.843, which is within 
the uncertainty of the measured ratio. 

It is apparent from Fig. QT] that for the QWL coating, 
two of the measured loss angles (the fourth and sixth) 
are somewhat high compared to the other six. For the 
optimized coating, one of the values (the eighth) is high 
relative to the others. Each of these three values is signifi- 
cantly more than one standard deviation above the mean 
for that coating. These anomalously high measurements 
may be due to transients that have been observed to sud- 
denly boost the noise floor of the instrument. They are 
not related to Brownian noise in the coatings. If these 
measurements are discarded, the means for the coatings 
are reduced slightly to the values indicated by asterisks 
in Table HI ^ QWL = 8.25 x 10~ 6 and opt = 6.85 x 10~ 6 . 
The ratio of the means then becomes r = 0.83 and the 
remaining loss angles for each coating are more tightly 
grouped. While there is no compelling justification for 
removing these points, the new means may reflect the 
true coating loss angles more accurately. 



It should be noted that there is an uncertainty asso- 
ciated with each value in Table [TT] since each is obtained 
from a least-squares fit to a spectrum with some scatter 
in it. This is evident in Fig. a plot of one of the 
spectra of the QWL coating overlaid by one of the op- 
timized coating. While there is some overlap between 
the spectra of the two coatings there is none between the 
confidence intervals of the fits (indicated by the dashed 
lines in the plot), because the error of the fit is much 
smaller than the standard deviation of the residuals to 
the fit. Figure [TBI is an histogram of the residuals of the 
fit to the QWL spectrum. The standard deviation of the 
histogram is 1.2 x 10 -19 m/\/Hz or about 10% of the 
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Dataset 

FIG. 11: Plot of the loss angles in Table El The black squares 
are the loss angles of the QWL coating, the black dots are the 
loss angles of the optimized coating. The solid gray lines are 
the means of each set of measurements, and the dashed gray 
lines are one standard deviation from the mean. 



length noise at 3 kHz and nearly equal to the separation 
between the spectra of the two coatings. The error of 
the fit, however, is smaller than the standard deviation 
by a factor of 1/ y/Nr>, where Njj is the number of data 
points. Finding the least-squares fit to the data is akin 
to finding the center of the histogram which can be done 
with much more precision than its standard deviation, 
particularly if there are many data points. In fact, this 
fitting error is smaller than the errors quoted above for 
the mean loss angles and can be ignored. 

One more thing to note from Fig. [13] is that the dis- 
tribution of residuals is consistent with the normal dis- 
tribution, suggesting that the scatter in each spectrum 
is random in nature and that our method of extracting 
the loss angles from least-squares fits to the length noise 
spectra is valid. 

The separation between the spectra of the two coat- 
ings in Fig. [T^] can be made visually more evident by 
using the Savitzky-Golay smoothing filter 

.15] 

to filter 

out the random noise in the spectra. This technique gen- 
erates a smoothed data set from the raw data where each 
point in the smoothed data is derived from a polynomial 
regression performed using neighboring points in the raw 
data. The smoothed data can be viewed as a generalized 
running average of the raw data. In the smoothed spec- 
tra shown in Fig. 1 141 each data point was derived from 
a first order polynomial regression using the correspond- 
ing point (at the same frequency) in the raw spectra of 
Fig. [12] plus the five points to the left and the five points 
to the right of that point. 



8 



QWL and Optimized Coatings 




2200 



2400 2600 2800 3000 3200 3400 36003800 
Frequency (Hz) 



FIG. 12: Plot of one of the spectra of the QWL coating (light) 
overlaid by one of the optimized coating (dark) . The solid line 
through each spectrum is the best fit of Eq. (10) to the data. 
For each fit, the region between the dashed lines represents 
the 95% confidence interval. 



CONCLUSIONS 



We have directly measured broadband Brownian noise 
of two different high reflectivity optical coating designs: a 
standard QWL coating, and an optimized coating specif- 
ically designed to minimize the Brownian noise. The 
ratio of the coating loss angles that we observed was 
r = 0.82 ± 0.04, which agrees with the predicted ratio, 
to within the margin of error. The results validate the 
proposed coating optimization strategy, and suggest its 
use to improve the sensitivity of future generations of 
interferometric gravitational wave detectors at relatively 
little cost. 
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FIG. 13: Histogram of the residuals of the fit to one of the 
spectra taken for the QWL coating. The distribution of the 
residuals is well approximated by a normal distribution with 
a standard deviation of 1.2 x 10~ 19 m/\/Hz. 
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FIG. 14: This is the plot in Fig. [12] after the Savitzky-Golay 
smoothing filter was applied to each spectrum. The separa- 
tion between the spectra of the two coatings can be seen more 
clearly. 



[1] P. Fritschel, Proceedings of SPIE 4856, 282 (2003). 

[2] Advanced LIGO Reference Design, 

LIGO document M060056, URL 

https : //dec . ligo . org/cgi-bin/DocDB/RetrieveFile?docid=1507 
[3] B. Willke et al., Class. Quantum Grav. 23, S207 (2006). 
[4] S. D. Penn et al., Classical and Quantum Gravity 20, 

2917 (2003). 
[5] R. Kubo, Rep. Prog. Phys. 29, 255 (1966). 
[6] G. M. Harry et al., Class. Quantum Grav. 24, 405 (2007). 
[7] Generated using BENCH now 

superseded by GWINC, URL 

http : //ilog ■ ligo-wa . caltech . edu : 7 285/advligo/GWIMC 
[8] K. Numata, Ph.D. thesis, University of Tokyo (2002). 
[9] E. D. Black et al., Phys. Lett. A 328, 1 (2004). 
[10] R. W. P. Drever et al., Applied Physics B 31, 97 (1983). 
[11] E. D. Black, Am. J. Phys. 69, 79 (2001). 
[12] S. J. Orfanidis, Electromagnetic waves and antennas, 

URL |http : / /www . ece . rutgers . edu/~ orf anidi/ ewaTj 
[13] G. M. Harry et al., Class. Quantum Grav. 19, 897 (2002). 
[14] J. Agresti et al., Proceedings of SPIE 6286, 628608 
(2006). 

[15] A. Savitzky and M. J. E. Golay, Analytical Chemistry 
36, 1627 (1964). 



